Internal states of model isotropic granular packings. II. Compression and pressure cycles

This is the second paper of a series of three investigating, by numerical means, the geometric and mechanical properties of spherical bead packings under isotropic stresses. We study the effects of varying the applied pressure P (from 1 or 10 kPa up to 100 MPa in the case of glass beads) on several types of configurations assembled by different procedures, as reported in the preceding paper. As functions of P, we monitor changes in solid fraction \Phi, coordination number z, proportion of rattlers (grains carrying no force) x0, the distribution of normal forces, the level of friction mobilization, and the distribution of near neighbor distances. Assuming the contact law does not involve material plasticity or damage, \Phi is found to vary very nearly reversibly with P in an isotropic compression cycle, but all other quantities, due to the frictional hysteresis of contact forces, change irreversibly. In particular, initial low P states with high coordination numbers lose many contacts in a compression cycle, and end up with values of z and x0 close to those of the most poorly coordinated initial configurations. Proportional load variations which do not entail notable configuration changes can therefore nevertheless significantly affect contact networks of granular packings in quasistatic conditions.Comment: Published in Physical Review E 12 page

Incremental response of granular materials: DEM results

We systematically investigate the incremental response of various equilibrium states of dense 2D model granular materials, along the biaxial compression path (\sigma 11 < \sigma 22, \sigma 12 = 0). Stress increments are applied in arbitrary directions in 3- dimensional stress space (\sigma 11, \sigma 22, \sigma 12). In states with stable contact networks we compute the stiffness matrix and the elastic moduli, and separate elastic and irreversible strains in the range in which the latter are homogeneous functions of degree one of stress increments. Without principal stress axis rotation, the response abides by elastoplasticity with a Mohr-Coulomb criterion and a non-associated flow rule. However a nonelastic shear strain is also observed for increments of \sigma 12, and shear and in-plane responses couple. This behavior correlates to the distribution of friction mobilization and sliding at contacts.Comment: 4 page

How granular materials deform in quasistatic conditions

Based on numerical simulations of quasistatic deformation of model granular materials, two rheological regimes are distinguished, according to whether macroscopic strains merely reflect microscopic material strains within the grains in their contact regions (type I strains), or result from instabilities and contact network rearrangements at the microscopic level (type II strains). We discuss the occurrence of regimes I and II in simulations of model materials made of disks (2D) or spheres (3D). The transition from regime I to regime II in monotonic tests such as triaxial compression is different from both the elastic limit and from the yield threshold. The distinction between both types of response is shown to be crucial for the sensitivity to contact-level mechanics, the relevant variables and scales to be considered in micromechanical approaches, the energy balance and the possible occurrence of macroscopic instabilitie

The nature of quasistatic deformation in granular materials

4 pagesInternational audienceStrain in granular materials in quasistatic conditions under varying stress originate in (I) contact deformation and (II) rearrangements of the contact network. Depending on sample history and applied load, either mechanism might dominate. One may thus define rheological regimes I and II accordingly. Their properties are presented and illustrated here with discrete numerical simulation results on sphere packings. Understanding the microscopic physical origin of strain enables one to clarify such issues as the existence of macroscopic elasticity, the approach to stress-strain relations in the large system limit and the sensitivity to noise

Flow of wet granular materials: a numerical study

We simulate dense assemblies of frictional spherical grains in steady shear flow under controlled normal stress $P$ in the presence of a small amount of an interstitial liquid, which gives rise to capillary menisci, assumed isolated (pendular regime), and to attractive forces. The system behavior depends on two dimensionless control parameters: inertial number $I$ and reduced pressure $P^*=aP/(\pi\Gamma)$, comparing confining forces $\sim a^2P$ to meniscus tensile strength $F_0=\pi\Gamma a$, for grains of diameter $a$ joined by menisci with surface tension $\Gamma$. We pay special attention to the quasi-static limit of slow flow and observe systematic, enduring strain localization in some of the cohesion-dominated ($P^*\sim 0.1$) systems. Homogeneous steady flows are characterized by the dependence of internal friction coefficient $\mu^*$ and solid fraction $\Phi$ on $I$ and $P^*$. We record fairly small but not negligible normal stress differences and the moderate sensitivity of the system to saturation within the pendular regime. Capillary forces have a significant effect on the macroscopic behavior of the system, up to $P^*$ values of several units. The concept of effective pressure may be used to predict an order of magnitude for the strong increase of $\mu^*$ as $P^*$ decreases but such a crude approach is unable to account for the complex structural changes induced by capillary cohesion. Likewise, the Mohr-Coulomb criterion for pressure-dependent critical states is, at best, an approximation valid within a restricted range of pressures, with $P^*\ge 1$. At small enough $P^*$, large clusters of interacting grains form in slow flows, in which liquid bonds survive shear strains of several units. This affects the anisotropies associated to different interactions, and the shape of function $\mu^*(I)$, which departs more slowly from its quasistatic limit than in cohesionless systems.Comment: 20 pages, 29 figures with 39 subfigure

Influences des paramètres micromécaniques dans la simulation numérique discrète des matériaux granulaires : assemblage, déformation quasi-statique, écoulements

27 pagesWe review the influence of micromechanical parameters on the macroscopic mechanical behaviour of granular materials, as numerically simulated in discrete element approaches, both in quasistatic conditions an din dense flow. We insist in particular on the role of suitably defined dimensionless numbers apt to provide a classification of rheological regimes of quite general validity

Quasistatic rheology and the origins of strain

Published in special issue of "Comptes-Rendus Physique" on granular materialsInternational audienceFeatures of rheological laws applied to solid-like granular materials are recalled and confronted to microscopic approaches via discrete numerical simulations. We give examples of model systems with very similar equilibrium stress transport properties -- the much-studied force chains and force distribution -- but qualitatively different strain responses to stress increments. Results on the stability of elastoplastic contact networks lead to the definition of two different rheological regimes, according to whether a macroscopic fragility property (propensity to rearrange under arbitrary small stress increments in the thermodynamic limit) applies. Possible consequences are discussed

Discrete numerical simulation, quasistatic deformation and the origins of strain in granular materials.

International audienceSystematic numerical simulations of model dense granular materials in monotonous, quasistatic deformation reveal the existence of two different régimes. In the first one, the macroscopic strains stem from the deformation of contacts. The motion can be calculated by purely static means, without inertia, stress controlled or strain rate controlled simulations yield identical smooth rheological curves for a same sample. In the second régime, strains are essentially due to instabilities of the contact network, the approach to the limits of large samples and of small strain rates is considerably slower and the material is more sensitive to perturbations. These results are discussed and related to experiments : measurements of elastic moduli with very small strain increments, and slow deformation (creep) under constant stress